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Abstract. Despite a Curie-Weiss temperature 9cw ~ —14 K, the Tb2Ti207 
pyrochlore magnetic material lacks long range magnetic order down to at least T* w 50 
mK. It has recently been proposed that the low temperature collective paramagnetic 
or spin liquid regime of this material may be akin to a spin ice state subject to 
both thermal and quantum fluctuations — a quantum spin ice (QSI) of sorts. Here 
we explore the effect of a magnetic field B along the [111] direction on the QSI 
state. To do so, we investigate the magnetic properties of a microscopic model of 
Tb2Ti207 in an independent tetrahedron approximation in a finite B along [111]. 
Such a model describes semi-quantitatively the collective paramagnetic regime where 
nontrivial spin correlations start to develop at the shortest lengthscale, that is over a 
single tetrahedron, but where no long range order is yet present. Our results show that 
a magnetization plateau develops at low temperatures as the system develops B = 
ferromagnetic spin- ice- like "two- in/two-out" correlations at the shortest lengthscale. 
From these results, we are led to propose that the observation of such a [111] 
magnetization plateau in Tb2Ti207 would provide compelling evidence for a QSI at 
B = in this material and help guide the development of a theory for the origin of its 
spin liquid state. 



1. Introduction 

The spin liquid (SL) is the state of a magnetic system characterized by short range spin- 
spin correlations down to temperatures much lower than the Curie- Weiss temperature, 
^CW; which is set by the magnetic interactions. While the SL concept was originally 
proposed in 1973 in the context of a spin-1/2 Heisenberg antiferromagnet on a 
triangular lattice, [1] it is only in recent years that its systematic experimental search 
in highly frustrated magnetic systems has really taken off. [2] A number of quasi two- 
dimensional [3] and three-dimensional (3-D) [H [5] materials with SL behavior have been 
tentatively identified. Yet, the Tb2Ti207 insulator is perhaps one of the 3D highly 
frustrated magnetic systems with SL behavior that has been attracting interest for the 
longest time. [5] In Tb2Ti207, the magnetic Tb'^"'' ions reside on a pyrochlore lattice 
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of corner-sharing tetrahedra. Despite ^cw ~ "14 K, [B] Tb2Ti207 does not develop 
conventional magnetic long range order down to at least T* = 50 mK. O [71 El E] 

It is unclear why Tb2Ti207 fails to order. [10] It was originally thought that the 
Tb'^"'" magnetic moments could be described by classical Ising spins that can only point 
in or out of a given tetrahedron, P, [H] similar to the Ho^"^ and Dy^+ moments in 
the Ho2Ti207 |T2j and Dy2Ti207 [13] spin ice materials. [H] However, such an Ising 
model for Tb2Ti207 predicts a transition to a four sublattice Neel ordered state at a 
temperature around 1 K, [TT] in dramatic disagreement with experiments. [3, [8] The 
microscopic justification for depicting the Tb^+ moments as Ising spins stems from 
considering only the single ion crystal field (CF) ground state doublet of Tb^+ and 
neglecting higher excited levels when describing the system at temperatures T < 10 
K. [6] However, a number of experimental [3, [HI [15] and theoretical [Ml [17] results argue 
against such an Ising description. For example, an Ising model fails to explain, even at 
a qualitative level, the symmetry of the neutron scattering intensity pattern, /(g), even 
in the paramagnetic phase. [3 HE] Interestingly, I{q) remains largely unaltered from 4 
K [7] down to 50 mK, never extending much beyond the size of a tetrahedron primitive 
basis. [8] One must therefore incorporate the excited CF states when constructing an 
effective low-energy Hamiltonian, i/efr, describing Tb2Ti207 at T < 10 K. The main 
reason for this necessity is that, unlike in spin ices, [18] the exchange and dipolar 
interactions in Tb2Ti207 are not much smaller than the energy gap A 18 K separating 
the ground and first excited CF doublets. [HI [3, [19] Also, when described as a classical 
(111) Ising model, Tb2Ti207 is found to be close to the boundary between an "all- 
in/all-out" four sublattice Neel ordered ground state and a spin ice state. [SI [11] Hence, 
perturbations, such as those introduced by virtual excitations to the excited crystal field 
states, [20] can in principle give rise to new ground states, and these excitations must 
therefore be considered at the very outset. 

A recent work [20] argued, on the basis of perturbation theory calculations that 
consider an independent tetrahedra approximation (ITA) to model the SL state of 
Tb2Ti207, that exchange and magnetic dipole-dipole interactions induce significant 
admixing between the ground and excited CF states. Most importantly, this admixing 
renormalizes the Ising (longitudinal) part of i/efr from that of a non-frustrated (111) 
Ising antiferromagnet [21] to that of a frustrated (111) Ising ferromagnet, in essence 
turning the system into a spin ice. [121 [3 [HI [221 [23] Interaction-induced admixing of 
the CF states also generates effective off-diagonal (transverse) spin-spin interactions in 
iJcfT- [201122] As a result, Tb2Ti207 may possibly be described by an effective frustrated 
(111) Ising ferromagnet model with additional quantum fluctuations transverse to the 
local (111) directions. [22] In other words, Tb2Ti207 may be viewed in the temperature 
range [50 mK < T < 10 K] as being in a thermally and quantum mechanically 
fluctuating spin ice state where frustrated ferromagnetic correlations in the (111) Ising 
components exist at the shortest lengthscale — a quantum spin ice (QSI) of sorts. [21] 
Hence, Tb2Ti207 may constitute an exciting playground to explore quantum effects in 
a spin ice system with the QSI proposal providing a long awaited and useful theoretical 
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perspective to rationalize Tb2Ti207. [TU] The possibility of a QSI in Tb2Ti207 is also of 
somewhat fundamental interest as it may relate to the topical problem of fractionalized 
spin excitations and emerging photon-like excitations in theoretical models of frustrated 
quantum Ising antiferromagnets on the pyrochlore lattice. [25] 

Because of the accidental quasi-cancellation of the exchange and dipolar energy 
scales at the nearest-neighbor level, the zero-frequency QSI correlations in Tb2Ti207 
may be masked by strong dynamical fluctuations and difficult to detect via neutron 
scattering [HI [TJ [HI E] and muon spin relaxation experiments. [5] From this perspective, 
the behaviors exhibited by Tb2Ti207 may be viewed as an interesting example of the 
spectral weight downshift expected in highly frustrated magnetic systems. [26] We 
believe that this phenomenology is an important part of the physics of Tb2Ti207 and at 
the origin of the stumbling blocks in developping a microscopic theoretical understanding 
of Tb2Ti207. [To] With this in mind, we explore here whether low temperature bulk 
magnetization, M{T), measurements may be used to expose whether the nearest- 
neighbor part of the longitudinal (Ising) sector of Hes describing Tb2Ti207 is indeed 
that of a frustrated (111) ferromagnet [21] and, therefore, whether a QSI description of 
this system is correct. Perhaps one of the clearest tell-tale bulk magnetic signatures 
of frustrated ferromagnetic "two-in/two-out" ice-rule correlations [H] in a spin ice 
material is the emergence of a magnetization plateau for a magnetic field B along 
the [111] direction (diagonal of the conventional cubic unit cell) at low temperatures 
once the ice-rule obeying spin correlations have become established. [271 128] Motivated 
by this magnetization plateau phenomenon in classical Ising spin ices, we therefore ask 
here whether a plateau-like feature in the [111] magnetization can also develop in the 
collective paramagnetic regime of Tb2Ti207 as the frustrated ferromagnetic (111) spin 
correlations develop at the shortest (single tetrahedron) lengthscale. 

The suggestion that Tb2Ti207 may be described at low energies and low 
temperatures by an effective Hamiltonian, Hes, whose longitudinal Ising components 
are frustrated by effective nearest-neighbor couplings was motivated by perturbation 
theory calculations carried out over a single tetrahedron of interacting Tb^"*" ions. [20] 
The coupling between the (111) Ising components for a pair of interacting Tb^"*" ions is 
in principle renormalized by the interactions with all the other Tb'^^ ions in the system. 
Therefore, an important question is whether the frustrated "two-in/two-out" spin-ice- 
like correlations are preserved in a perturbation theory extended to a thermodynamically 
large lattice and, in particular, in the presence of the long-range magnetostatic dipole- 
dipole interactions between Tb'^+ ions. [6l [TTl [16] Our preliminary perturbation theory 
calculations carried over the whole lattice show that the system remains in an ice-rule 
obeying ground state with competing states with a "three- in/one-out" spin configuration 
having a slightly higher energy. [20l [22] Consequently, we would expect that, in an 
"exact" treatment of Hes at nonzero temperature, spin ice "two-in/two-out" correlations 
would develop at low temperatures. The important question is therefore whether such 
correlations actually develop in real Tb2Ti207 at low temperatures and if they can be 
evinced by the observation of a [111] magnetization plateau. 
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To get a preliminary handle on the relevant temperature and magnetic field scale 
where QSI physics may become experimentally discernable in Tb2Ti207, we consider 
the same ITA used in Ref. [20] to describe the SL state of that material since, as stated 
above, interactions beyond a single tetrahedron do not spoil the existence of an ice-rule 
obeying semi-classical ground state. [22] The merit of the ITA to describe the collective 
paramagnetic state of Tb2Ti207 was discussed in Refs. [201 [221 [29] . In the context of the 
work presented here, one could interpret the ITA as a self-consistent cluster mean-field 
theory (SCMFT) [30] where the considered cluster is a tetrahedron, hut where the self- 
consistency conditions that incorporate the inter-clusters (inter-tetrahedra) interactions 
have been set to zero. This is because we are describing the collective paramagnetic/spin 
liquid regime of the model where the correlations extend solely over the smallest cluster 
lengthscale where no global symmetry is (yet) spontaneously broken. The experimental 
observation that spin correlations in Tb2Ti207 never develop beyond the size of a 
single tetrahedron from a temperature of ~ 10° K down to 50 mK [5l [71 [H [9] would 
seem to further provide a post-factum justification to this "constrained" paramagnetic 
description of a SCMFT of the H^k of Tb2Ti207. [31] In other words, the ITA should 
be qualitatively valid for sufficiently small field and temperatures such that long-range 
spin-spin correlations are not induced. [15] Finally, for the proposed experiment seeking 
to expose a magnetization plateau for B along [111], symmetry considerations play a 
crucial role, [271 [2H] and the ITA does capture that symmetry aspect of the problem. 
Below, we find that for parameter values for which the model exhibits a QSI, and which 
are appropriate for Tb2Ti207, [20] the low-temperature magnetization M(T) (T < 100 
mK) does indeed exhibit a plateau for B = \B\ ~ 0.1 T. 

2. Model and Results 

Within the ITA, the Hamiltonian for Tb2Ti207 reads 

a=l,b>a a=l 

Hf^ = I)/2i[J„- Jfe-3(J,-f,6)(Jfe-f,,)]|i?,,|-3 and = JJa ■ Jt describe, 

respectively, the dipole-dipole and the exchange interactions for a pair ab. = 
—gUsJa ■ B is the Zeeman Hamiltonian and H^^ is the crystal field (CF) Hamiltonian, 
both for ion a. [HI [19] Tiab = Rfe — R-a = iR-abl'^afe where Ha is the position vector of ion 
a, and -Rnn = 3.59A is the nearest-neighbor distance. V = /io(5'/iB)^/(4vri?nn^) = 0.0315 
K is the dipole-dipole coupling and J' = 0.167 K is the exchange coupling, with the 
convention that ^7 > is antiferromagnetic. [SO, [32] /xb is the Bohr magneton, g = ^ is 
the Lande factor of Tb'^+ and (| J^l =6) is the total angular momentum operator of 
Tb^"*" ion a. Calculations considering a single (non-interacting) Tb^+ ion in nonzero B 
show that the six lowest energy CF states describe well the exact M when compared 
with calculations that consider all 2J+1=13 CF states for i? < 30 T with an error 
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less than 1%. Henceforth, we use the six lowest CF states of Tb^"*" to calculate M 
within the ITA. The basis states of the system are taken as the tensor product of 
the CF states of the four non-interacting Tb^"*" ions on a tetrahedron. We denote the 
two states of the CF ground state by | ]) and | [). These correspond, respectively, 
to the "out" and "in" degenerate (Ising) CF ground states of = Yla 
"upward pointing" tetrahedron primitive basis [6], [11], [27j when referring to the [111] 
direction along which B points. The CF ground doublet and first excited doublet 
are given in Ref. [20] and the two other excited states, determined as in Ref. [20], are 
= a6(|6) + |-6)) + a3(|3)-|-3))and l^^""^) = be{\Q) + \-6)) + h{\3) -\-3)) + bo\0) 
with energies E-s = 142.1 K and E4 = 212.4 K, respectively. [33] Here we express \ip^^'') 
and \tp^^^) in terms of the eigenstates | J = 6, mj) of within the fixed J = 6 manifold. 
The local z easy axes for sublattices 1 to 4 on an individual tetrahedron primitive unit 
cell are along [111], [111], [111] and [111] and we maintain this sequence of labelled 
spins/sites throughout the paper, with \ijkl) = \ijkl)i234 = ® \j)2 ® Ik)^ ® |/)4. 

We determine the matrix elements of the full H in Eq. (1) in the CF basis, find 
the eigenvalues and eigenstates of H using a numerically exact diagonalization method 
and use these to calculate M given by the standard formula 

Ns 4 

= 3J^J^J2<'iM\<P^)nBiE,) , (2) 

i=l a=l 

where M° is the a component of the total M per ion. |0j)'s are the eigenstates 
of Hamiltonian ([T]), uf is the a(3 component of the rotation matrix from the local 
quantization frame, aligned along the local cubic [111] direction at site a, to the global 
axis. J° is the a component of J expressed in the local xyz quantization frame with z 
aligned with the local [111] axis. n-B{Ei)=exp{—Ei/kBT) is the Boltzmann population 
of state i with energy Ei at temperature T. The sum over i in Eq. (2) is carried over 
the A'^s = 6^ = 1296 eigenstates retained in H'^^ and the sum on a is over the four sites 
of a tetrahedron. 

We calculate M within the ITA for B along [111] and for T=20, 50 and 100 mK. 
The result of this calculation is shown in Fig. [1] As in spin ices, [27] a magnetization 
plateau occurs, becoming apparent only below the temperature at which the ice rules 
are well established. [20] That is, when the temperature drops below the very small 
energy gap separating the manifold of "two-in/two-out" states from the manifold of 
the "three- in/one-out" states. In the present model, the spin-ice phenomenology is 
"hidden" at low temperature because of the accidental location of the system near the 
phase boundary between the "all-in/all-out" (doublet) Neel state and the "two-in/two- 
out" spin ice state. [20] The plateau develops only at low temperature because the QSI 
state in H in Eq. (1), and possibly in Tb2Ti207, is inherited by perturbative virtual 
quantum mechanical crystal field excitations renormalizing the already accidentally 
almost vanishing classical (longitudinal Ising) interactions in Hgg. [20] Hence, the 
manifold of low-energy states at B = 0, which descend from the 2^ = 16 parent Ising 
CF ground doublet states, span a small energy bandwidth 6Q ~ 0.5 K, [20] with the 
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B(T) 

Figure 1. (Color online). Magnetization M of Tb2Ti207 in the ITA as a function 
of magnetic field B applied along the [111] direction. In the magnetization plateau 
regime, where M{B) flattens out for T = 20 mK, the system is in a "quantum kagome 
ice" state (see text). Inset: energy of the lowest states as a function of B which which 
shows a level-crossing at i3c=0.082 T. 

6 lowest energy spin ice states spanning an even smaller bandwidth 6uj ~ 0.06 K (see 
the six lowest [singlet © doublet © triplet] states at -B = in the inset of Fig. 1). [20] 
Consequently, the plateau in M only emerges at quite low temperature. This constitutes 
a prime example where high frustration is caused by naturally accidentally fine-tuned 
microscopic interactions (here H'^^ and H'^^^) and where the relevant low-energy physics 
is pushed down to an exceedingly small scale. [26] For small B, M starts to increase 
even at a very low T (e.g. 20 mK). The reason for this is that, despite the ground state 
being a singlet at the ITA level, there exist a second order (van Vleck like) susceptibility 
originating from the excited states which leads to M 7^ for B > 0. 




B(T) 

Figure 2. (Color online). [Ill] magnetization M of Tb2Ti207 in the independent 
tetrahedron approximation as a function of B in the [111] direction. Here the 
temperature is 20 mK. 

The value of the exchange J' in Tb2Ti207 is not fully agreed upon. [51 [32] It is 
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therefore worthwhile to explore how M{B) changes as is varied. We plot in Fig. [2] 
M as a function of B for T= 20 mK and exchange couplings = 0.15, 0.167, 0.18, 
0.19, 0.20 K. These results show that a magnetization plateau is manifest only when 
the system displays a QSI state at B = ioi J' < J^. [20] Here J7c=0.187 K delineates 
the transition between the QSI and the "all-in/all-out" (doublet) state, [20] with the 
plateau disappearing for J > J^. For J = 0.18 K, the system is barely beyond the 
threshold for exhibiting a QSI state and one finds that the plateau is poorly formed 
even at T = 20 mK. Our calculations also show that, although a rescaling of and V 
changes the magnetization plateau value and the B range over which a magnetization 
plateau arises, it does not change the generic qualitative behavior seen in Fig. 1 as long 
as the model, and presumably the real system, displays a QSI state for B = or, in 
other words, as long as the system possesses frustrated ferromagnetic interactions in the 
nearest-neighbor Ising part of Hes- 

3. Field Dependence of Energy Levels on a Single Tetrahedron 

The ITA results above can be interpreted as describing the behavior of the system in the 
temperature range where no long range order exists except for the field-driven global 
[111] magnetization, but where spin correlations extend over the shortest lengthscale 
characteristic of a collective paramagnetic/spin liquid regime (i.e. over a tetrahedron). 
Yet, it is interesting, for completeness, to consider the evolution of the energy levels over 
a single tetrahedron as this evolution would ultimately need to be taken into account in 
a SCMFT [30] going beyond a single tetrahedron to describe the possible spontaneously 
broken symmetry states of the model. 

To understand M{B) at 20 mK within the ITA, we scrutinize the ground state 
of the system as a function of B. We first note that the ground state, |\E'), in the 
B range considered in Fig. [1], can be generically and compactly written as |^) = 
a\u)+a\u)+(3\nU) + (^\x)- Here \uj) = \ ttii) + | Titi) + I tiit) and \u) is the 
time conjugate of lo;). a, a, P, e and |x) are 5-dependent and e is much smaller than at 
least one of the coefficients a, a, or /5. Notice that in la;), ion 7^1 is in the 1 1) ("out") 
state. We denote the QSI state with B = in Eq. (1) as |QSI) = l^') above with a = a, 
P = and |e| <^ 

The evolution of the ground state with B can be explained from the competition 
between and H^'^' + H'^'p in Eq. (1). For weak B, {QSl\H^\(^) ^ 0, where \(p) refers to 
any of the six pure "two-in/two-out" states (e.g. Ittii)) I Titi)) etc). Among these, the 
three states with ion 7^ 1 in| Di (i.e. polarized along [111]) have a lower Zeeman energy 
than the other three two-in/two-out states. Hence, with B ^ 0, a state with larger 
weight from these components, i.e. with ion #1 | has lower energy than the B = 
I QSI) state and become the ground state. As a result, for non-zero and sufficiently weak 
B, the ground state evolves from |QSI) to a state |\E') with a > a and with = (/3 = 
because does not admix |QSI) with | tiii) at any order in H^). 

As B is increased, a grows and a decreases until B* ^ 0.016 T where dE"^ = 
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(QSI|i/^|QSI) ~ 0.06 K. For B ~ i?*, all 6 lowest energy states with a predominant 
two-in/two-out weight spanning the small energy bandwidth 5uj ~ 0.06 K (see inset 
of Fig. 1) fall within the energy range 6E'^ . Hence, for B > B*, the system adopts a 
ground state with a ^ a, where ion # 1 is polarized towards [111], while the three other 
ions are in an entangled state which is a symmetric linear combination of two- in/one-out 
states. This partially spin polarized (PSP) state, |PSP), is analogous to the kagome ice 
state in classical spin ice systems p7l [28] . By increasing B beyond B*, up to B^, the 
ground state remains |PSP) (i.e. P = because does not admix |PSP) with | till) 
at any order in H^), until a level crossing occurs for B^ ~ 0.082 T (see inset of Fig. [1]). 

One must consider the evolution of the first excited energy level with B in order 
to understand the above level crossing. For B = 0, one of the two states from the 
first excited doublet has a small | tiii) "three-in/one-out" component. As a result, the 
matrix elements of between this state and the higher energy excited states which 
have a three- in/one-out component is nonzero. Hence, by increasing B, the | Tiii) 
weight in the first excited state increases and its energy decreases and, by the time 
B ~ 0.08 T, this state is almost pure | tiii)- This state then crosses the PSP state 
at Be ~ 0.082 T and becomes the ground state (see inset of Fig. ([T]), with the system 
exiting the PSP state characterized by a magnetization plateau. 

4. Conclusion 

We calculated the magnetization M of Tb2Ti207 for a field B along the [111] direction 
within a single tetrahedron approximation (ITA). For weak B in this direction, a 
magnetization plateau is found to occur for a temperature range where the system 
is in a collective paramagnetic/spin liquid regime. Within this regime, the correlations 
extend solely over the size of a tetrahedron and each tetrahedron primitive basis is 
in a state predominantly weighted by |t)i ® (|Tii)234 + |iTi)234 + |iiT)234), where the 
three spins on the vertices of the triangle perpendicular to the [111] direction, i.e in 
the kagome plane, form a linear combination of two-in/one-out states. [271 128] This 
partially spin-polarized state (PSP) is reminiscent of the classical kagome ice state in 
classical Ising spin ices [271 128] which one may label quantum kagome ice to distinguish 
it from the degenerate kagome ice state. From a microscopic level point of view, both 
the zero field quantum spin ice (QSI) and the PSP states originate from virtual crystal 
field excitations and quantum many-body effects. [20] At this time, low-temperature 
experiments are needed to test the proposals of QSI and PSP states in Tb2Ti207. From 
the results presented here, we propose that an experiment based on measurements of 
M along the [111] direction would help to ascertain the possibility of a QSI state in 
Tb2Ti202. Furthermore, the observation of a QSI and a PSP in Tb2Ti207 would confirm 
the status of this material as a unique opportunity to explore quantum effects in a spin 
ice system. It may also open an experimental avenue to explore some of the exciting 
physics proposed to be at play in frustrated pyrochlore Ising systems with perturbative 
quantum fluctuations. [25] We hope that this work will motivate such investigations 



Proposal for a [111] Magnetization Plateau in the Spin Liquid State of Th2Ti20'j 9 

which would help shed some light on the microscopic origin of the enigmatic spin liquid 
state in Tb2Ti207. [3 [10] 
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